Entanglement in an SU(n) Valence-Bond-Solid State 
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We investigate entanglement properties in the ground state of the open/periodic SU(n) 
generalized valence-bond-solid state consisting of representations of SU(7i). We obtain ex- 
act expression for the reduced density matrix of a block of contiguous spins and explicitly 
evaluate the von Neumann and the Renyi entropies. We discover that the Renyi entropy is 
independent of the parameter a in the limit of large block sizes and its value 2 log n coincides 
with that of von Neumann entropy. We also find the direct relation between the reduced 
density matrix of the subsystem and edge states for the corresponding open boundary system. 
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INTRODUCTION 



There is considerable current interest in quantifying entanglement in various quantum many- 



body systems. Entanglement in spin chains, corre! 
models was studied in detail [l|, U, \i, U, U U U, M, 



ated e 
9, 10 



ectrons, interacting bosons and other 
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la, llSl: 



Entanglement is a resource for quantum computation, it shows how much correlation we can use 
to control quantum devices [2l|, |22]. There are several different measures of entanglement. The 



most famous is the von Neumann (entanglement) entropy of a subsystem 2]|. This measure has 
recently been used to detect quantum phase transitions and topological/ quantum order [23|, |2j] in 
strongly correlated systems. We can also use the Renyi entropy to quantify the entanglement. The 
Renyi entropy was first proposed in information theory 25|. The von Neumann entropy S{pa) and 
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the Renyi entropy Sa{pA) are defined as follows: 



S{pa) 

SaiPA) 



-Tr(pAlogpA) 
^ ■logTr(pl) 



a 7^ 1 and a > 0. 



(1) 
(2) 



1-Q 

Here pA is the reduced density matrix of subsystem A and the power a is an arbitrary parameter. 
The Renyi entropy characterizes the mixed state much better: if we know the Renyi entropy at 
any a we know all eigenvalues of the density matrix. 

Studying entanglement also helps to understand the physics of quantum spin systems [l|, y] . 



M 



23] plays a very 



The model introduced by Affleck, Kennedy, Lieb, and Tasaki (AKLT model) 

important role in condensed matter physics. The exa ct g round state of this model is known as 



the Valence-Bond Solid (VBS) state. In 1983, Haldane 28(] conjectured that the antiferromagnetic 
Heisenberg Hamiltonian describing half-odd- integer spins is gap-less, but for integer spins it has 
a gap. The AKLT model agrees with this conjecture and enables us to understand the ground- 
state properties of gapped spin chains in a unified fashion. The construction of the AKLT-type 
model is not restricted to one dimension. In fact, the AKLT model was formulated on an arbitrary 
graph and an integration over classical spins was used for the evaluation of correlation functions 
in the VBS ground state [29|. The VBS state has attracted revived interest from the viewpoint 
of quantum information theory. An implementation of the AKLT model in optical lattices was 
proposed recently [30|, and the use of the AKLT model for universal quantum computation was 



32( 1 and to 



discussed in |31[. The VBS state is also closely related to the Laughlin wave function 
the fractional quantum Hall effect [33|]. Entanglement in S* = 1 AKLT model was first considered 

3 511. An interesting 



in [Sj]- Then the results were generalized to the arbitrary integer spin case in 
generalization of the AKLT model to the SU(n) version was constructed in [36l . 

In this paper we study entanglement in SU(n) version of the AKLT model. We consider entan- 
glement of a block of spins with the rest of the ground state. We evaluate the von Neumann and 
the Renyi entropies of the block. We first confirm that the von Neumann entropy of a large block 
of spins reaches saturation. This is a partial proof of the conjecture proposed by Vidal et al [3]. In 
the SU(2) AKLT and the XY spin chains, this conjecture has alr eady been proved and the limiting 



n 



34|, 



35|, 



3d, l39|, y]. We find the 



entropy of the large block of spins was explicitly calculated a, 
essential simplification in the limit of large block sizes. In this limit, we discovered that the Renyi 
entropy is independent of a, actually it coincides with the value of von Neumann entropy 2 Inn. 
This means that the density matrix of the block is proportional to identical matrix of dimension 
n^. It is much different form XY spin chain, where the density matrix of the large block is infinite 



dimensional and eigenvalues are different (see 



4ll]). We also explore the connection between the 



reduced density matrix of the subsystem and degenerate ground states for the corresponding open 
boundary system. 

The paper is organized as follows. In the next section, we will study entanglement in the ground 
state of the SU(n) AKLT model with boundary spins. This section is the main part of this paper. 
The von Neumann entropy and the Renyi entropy will be evaluated for the block of neighboring 
spins. We will also discuss the direct relation between the reduced density matrix for the block 
and the degenerate ground states of the AKLT model with an open boundary condition. In the 
third section, we will investigate entanglement in the ground state of the SU(n) AKLT model with 
a periodic boundary condition. In this section, we will obtain the explicit form of the reduced 
density matrix. The finite size effect of the von Neumann and the Renyi entropies will be studied. 
The last section will be devoted to summary and discussions. 

II. SU(n) VBS STATE WITH BOUNDARY SPINS 
A. Construction of SU(n) VBS state 

In this section, we consider the SU(n) VBS state with boundary spins and calculate the von 
Neumann entropy (entanglement entropy) and the Renyi entropy of a block of contiguous spins. 
What we mean by "spin" in our system is an adjoint representation of SU(n). We shall first 
construct an SU(n) VBS state which consists of N adjoint representation of SU(n) in the bulk 
and fundamental and conjugate representations of SU(n) on the boundary. First, we prepare 
sites k {k = 0,1,..., A^) and k {k = 1,2,...,N + 1) and arrange SU(n) singlets consisting of a 



fundamental (D) and its conjugate (D) representations as shown in Fig. ([TJa) [43] ■ We assign \j) 
(j = 0, 1,2, ...,n — 1) to the fundamental representation, while \j) {j = 0, l,2,...,n — 1) to the 
conjugate representation. \j) can be represented by the tensor product of (n — 1) |j)s as 

|j)^-=i== Yl e^"^'-'"1«2,. ..,«„) (3) 

V(^-l)!a2,...,a„ 

where gJ<^2, ■■■,"« jg g, totally antisymmetric tensor of rank n. Using \j) and \j), an SU(n) singlet 
state l^o.o) can be represented as a maximally entangled state: 

N-l 

The above relation can be easily confirmed by inserting the resolution of the identity 1 = 
S?=o |j)OI s-iid substituting Eq.([3]). Next, we prepare the adjoint representation of SU(n) by pro- 



a) 
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b) 



(fundamental) 



n-1 



(conjugate) 



e 



}n 



(singlet) 



(adjoint) 
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FIG. 1: a) Construction of the SU(n) VBS solid state. White and black dots represent the SU(n) fundamental 
and its conjugate representations, respectively. A dimer corresponds to the singlet state I'I'co) and a circle 
denotes the projection onto the adjoint representation, b) The decomposition rule for the tensor product of 
D (g) D. c)Young tableaux corresponding to [n, n] and [n, n — 1, 1] 

jecting the tensor product D^D onto an (n^ — l)-dimensional subspace. This procedure corresponds 
to circles in Fig. [Ha). In Fig. [Hb), we visualize the decomposition rule Di^D = (singlet) ©(adjoint). 
Then we obtain the SU(n) adjoint representation at each composite site {k, k). Henceforth we shall 
call this composite site k. Finally, we can represent the SU(n) generalized VBS state as 



\G) 



>.N 



Pi.-,)\^ 



k=l^kk)\^0,0/0l\^0,0/12 ■ 



^r 



1$ 



0,0/ArAr+i5 



(5) 



where P/^-^ is a projection operator onto an adjoint representation of SU(n). 

Let us now consider the Hamiltonian whose ground state is the above state we have constructed. 
The Hamiltonian can be constructed along the same line as the AKLT model: 

N-l 

H =^Hk + 7ro,i + TTAf^AT+i, Hk = ^ ayPk^k+i^ (6) 

fc=l Y 

where y is a Young tableau which is neither [n, n] nor [n, n— 1, 1]. Here we have assigned [ki, ..., k\^\ 
to the Young tableau y, where Kj is the number of boxes in the j'-th column and Ai is the number 
of boxes in the first row. -P^^_)_x is a projection operator which projects (adjoint)(8)(adjoint) onto a 
representation characterized by Y and the coefficient ay can be an arbitrary positive number. The 
reason why [n, n] and [n, n — 1, 1] are excluded from the sum is following. Since D at site k and D at 



site fc+1 have already formed a singlet in the ground state ([5]), the possible representations obtained 
from the decomposition of (adjoint)® (adjoint) are restricted to [n,n] and [n, n — 1, 1] (graphically 
shown in Fig.l.c)). vro^i and tt^^n-^i are boundary terms which assure the uniqueness of the ground 
state of this Hamiltonian. ttq^i and vtjv^tv+i can be written in terms of the projection operators 
acting on the tensor products (fundamental) (^ (adjoint) and (conjugate)®(adjoint), respectively. 
By construction, the SU(n) VBS state ([5]) is a zero-energy ground state of this Hamiltonian. We 
should note here that another construction of Hamiltonian by Greiter and Rachel 



371 ] is similar but 



slightly different from ours. 



B. Reduced density matrix and von Neumann and Renyi entropies 

Next, we consider the reduced density matrices of subsystems of the ground state \G). To 
calculate the reduced density matrix, it is more convenient to recast the chain of singlets 
|^o,o)oil^o,o)i2 • • • l^o,o)jv7v+T ™ Eq.dni) in a different form. Let us first consider a chain of two 
singlets |<&o,o)oil*^o,o)i2- We can rewrite this product state as 

n— 1 n—1 

l^0,0>0ll^0,0>12 = -J2Yl l^',"^>02l^/,-m>Il, (T) 

1=0 m=0 

where \^m,n) is a basis of the maximally entangled state defined by 

Here / is an n-dimensional identity matrix and Ui^m = X^Z"^ {m, n = 0, 1, ..., n — 1) are generalized 
Pauli matrices, where the unitary operators X and Z act on \j) as X\j) = \j + l(mod.n)) and 
Z\j) = uj^\j) with u = ^'^'^1'^ ^ respectively. One can easily show the relation ([7]) by using the fact 
that |<l*o,o) is invariant under the action of iJJi^m ® ^i-'m) [43]. This procedure can be regarded as 
a multi-dimensional generalization of entanglement swapping. Repeatedly using the relation ([7|), 
we can generalize in a straightforward way to a chain of singlet states: 

|^0,o)oil^O,0>12--- I^O,o)Ar7v+T = ;^ X] ■■■ X] l^'i -"^l)Il••• l^«iv,-mjv)iVAr 

(il,mi) (ijv.miv) 

x(C//i,mi • • • ^l^,m^ ® -f)|$0,o)olV+T' (9) 

where {mkyrik) {k = 1,2,...,N) runs from (0,0) to (n — l,n — 1). To obtain the ground state 
\G) from ([9]), we have to make a projection onto the subspace of adjoint representation at each 
site k. Since the decomposition rule 0(8)0 = (singlet) © (adjoint) and the fact that |*l'o,o) is an 
SU(n) singlet, the vector space of the adjoint representation is spanned by \^i^-„b) ((^w-) / (0)0)). 



Then the only thing to do is to omit the summation over {Ik^fnu) = (0,0) in Eq.Q. The SU(n) 
generahzed VBS state can be rewritten as: 

1^*) = ,2 _ i\N/2 zZ '" zZ l^h,-mi)ll--- |^«]v-miv>iVAr(f^h,mi •••f^iiv.mjv ®-?')|^0,o)olv+T' 

7^(0,0) #(0,0) 

(10) 

where we have aheady normahzed \G) by the factor l/(n^ — 1) '^. 

Let us now consider the reduced density matrix of a block of contiguous spins of length L. 
We first suppose that the block starts from site k and stretches up to fc + L — 1, where k > 1 
and k + L — 1 < N. The reduced density matrix is obtained by taking the trace over the sites 



j = 0, 1, ..., k-l and j = k + L, ..., iV, iV + 1 as 

PL = ^^l,...,k-l,k+L,...,N,0,N+l\^)(^\ 

^ (^2 _ i\N 2^ " ' 2^ A^ " ' 1^ 1^ A^ ■■■ 2^ 2^ 

(Il.mi) (!fc_i,mj._i) (lf.^i^,'mf._f_^) {ijv™jv) ('fc.™fc) ('i>'"',,) ('L + fe-1 .'"i + fe-l) Cr , t._i >™'r , i._i ) 

7^(0.0) 7^(0,0) 7^(0,0) 7^(0,0) 7^(0,0) /(o,o) 7^(0,0) '-+" =^(0,0) 



L+fc-1' '"i+fe-l' 



X l^/fc -mjfcfc(% -m'J • • • \'^lL+k-l-mL+k-l)L+k-lL+k-l('^l 

X Tr^^jr^{U,VU2CSl)\'^ofi)o-w+i{^o,o\{UiV'U2CSl)^, (11) 

where Ui = Ui.^rm ■ ■ ■Ui^_^^mk_^, U2 = Ui^^^^m^^^^ ■ ■ -Ui^^m^, V = Ui^,m^ ■ ■ ■Ui^^^_^^mL+,_i and 
^' = ^i' m' ■ ■ ■ ^i' m' ■ To rewrite Eq. (fTTI) . we use the following property of l^o 0): 

fc'fc L + fc — l'Z/ + fc — 1 ' 

(5 r)|$o,o) = {ST' /)|$o,o) = (/ ® TS')\^o,o), (12) 

where S* and T are n-dimensional unitary operations acting on \j) and |j), respectively, and the 
superscript t denotes the transposition. Using this property and the cyclic property of the trace, 
we can simplify the last part of Eq. (fTTj) as 

^o,NTT(Ui^U2 ® I)\^o,o)ow+T{^o,omiV'U2 I)^ 
= Ti^j^iUi I){V I){U2 /)|$o,o)oivTT(^o,o|(f/2 ^ I)HV' I)HUi /)t 
= Tr^^jr+iiV /)(/ t/^)|<J>o,o)olvTT(^o,ol(/ 55 [/^)^(y' /)^ 
= Tro_j^(/ C/*)(y /)|<I>o,o)o1v+t(^o,oI(^' ® /)^(/ ® f^l)^ 
= '^o,N+iiy^I)\^ofi)oN+i{'^ofi\{V'(S>I)^. (13) 

Since ([13]) does not depend on {li,mi),- ■ ■ , (^fc-i, "^fc-i) and (4+L,'"T'fc+L), " " " ,(^Af,^Af), we can 



rewrite Eq. (jlip as 

7^(0.0) 7^(0,0) 7^(0,0) ^ #(0,0) 

From the form of the reduced density matrix p4p , we immediately notice that the reduced density 
matrix does not depend on both the starting site k and the total length of the chain A^. The same 



property for SU(2) 5 = 1 VBS state has already been proved in Ref. 3J|. We can regard above 
result as an SU(n) generalization of their result. 

Since the reduced density matrix is independent of both k and A^, we can set N = L without 
loss of generality. We can further reduce the original problem to that of the reduced density matrix 
of two end spins (D and D) using the following property of a bipartite pure state. Suppose that 
|^)a_b is a bipartite pure state of a total system AB. Then there exist orthonormal states \iI^j)a 
for the subsystem A, and orthonormal states \4'j)b for B such that 

|^>AB = Y^ ^Mj)A\(t>j)B. (15) 

3 

where Pj{> 0) satisfy YljPj ~ ^- This decomposition is called the Schmidt decomposition. The 



proof of the above theorem using the singular value decomposition can be found in Ref. 4j]. Prom 
Eq. (fT5]) . one can immediately notice that the set of eigenvalues of pA = Tii:b\'^)ab{'^\ coincides 
with that of ps = Tr^|^)^B(^|- 

Now we can reduce the eigenvalue-problem of pi to that of the reduced density matrix for end 
two spins pj^. p^ has the following form: 

PL = („2 ]_ i)L E ••• E (f/®/)|$o,o)o,X+T(^o,o|(C/®/)n (16) 

7^(0,0) 7^(0,0) 

where U = f/ii,mi ■ ■ ■ Uij^^rn^- To evaluate the eigenvalues of p^, it is convenient to formulate the 
action of {Ui^m tX" /) as a transfer matrix. Let us first see the action of {Ui'^m' ® -^) on a state \^i^m)'- 

(U,,^, I)\<^l,m) = (X''Z™'X'Z"^ /)|$o,0> = i0"'''\^l+,,m+n.'), (17) 

where we have used the relation ZX = ujXZ. Using the above relation, we can prove that 

iUl',m' f^ I)\'^l,m){'^l,m\iUl',m' (E) I)^ = \'^l+l' ,m+m') {'^l+l' ,m+m'\- (18) 

Next, we assign the vector (0, ...,0, l((/,m)-th entry), 0, ..., 0)* to the state \^i,m){^i,m\- This one 
to one correspondence plays an essential role in our analysis. From this bijection, the operation 



^ni m')jt(oo)(^i',rn' "SD I)\^i,m) {^i,m\{Ui' ^m' 'i^ I) ^aii be Written in terms of (n x n )-dimensional 
matrix as 



n 



1 1 

1 1 
1 1 



• • 1 

• • 1 

• • 1 



(19) 



\i 1 1 ... oy 

This transfer matrix can be diagonalized by the following unitary matrix: 



n 



(l 1 

1 c 

1 e 



1 

e 



1 /""^"l /-2(n2_i) 



1 



\ 



*2(n2-l) 



^{n2_i)2 y 



(20) 



where C = exp(27ri/n^). Then we can obtain the explicit form of the reduced density matrix p^ as 



PL 



-T^(1,0,...,0)* 



-C/,[diag(n2 - 1, -1, .., -l)fU^Ah 0, ..., 0)* 



(„2_l)L V >"'•••>"/ (n2-l)i 

= ;^(1 + ("'-1KW)I^o,o)o1vTt(^o,o| + ^ Y1 (l-Pn(^))l^',-)olvTT(^',™l'(21) 

(/,m)7^(0,0) 

where we have used the relation, ! + ('' + (^^ + . . . + C*-"" ~^^^ = (1 < /c < n^ — 1) and Pn{L) = 
{-TZi)^ ■ Substituting n = 2 into Eq. ([2T]) . one can reproduce the result of the SU(2) S=l VBS 



state obtained in Ref. 



m- 



Let us now start the evaluation of the von Neumann and the Renyi entropies of a block of L con- 
tiguous spins. First, we shall examine the von Neumann entropy of the block. From the Schmidt de- 
composition and the definition of the von Neumann entropy S{pl) = S{p£) = —Ty:i,2,...,LiPi I^SPl)' 
we obtain 

l + (n2-l)p„(L)^ 



5(pL) = 21ogn- 



n' 



-log(l + (n^-l)p„(L))-(n^-l) 



l-Pn{L) 



n^ 



log(l-p„(L)) (22) 



with Pn{L) = ij^T^) ■ Similarly to t 



chains in the gapped regime 



6 



34 



35| and the XY spin 



le SU(2) integer-5 VBS states 
3^, l39|, |40|] , S{pl) is bounded by 21ogn in the limit of large block 



sizes (L -^ oo) and approaches to this value exponentially fast in L. This is a partial proof of the 
conjecture proposed by Vidal et al. [3], that the von Neumann entropy of a large block of spins 



in gapped spin chains shows saturation. Next we shall examine the Renyi entropy of our system. 
From the definition of the Renyi entropy Sa{pL) = jr^ logTr(p2) (c^ 7^ l,and a > 0), 

SaiPL) = Y^ log(Ao,o(^)" + (n^ - l)A,,„^o,o(^)"), (23) 

where 

X .^^ i ^{l + {n^-^)Pn{L)), {l,m) = (0,0) 
[ J^(l-p„(L)) (/,m)/(0,0). 

Now we consider the limit of large block sizes, i.e., L — > 00. In this case, Xi^m become degenerate 
and great simplification occurs: 

SaiPa) = 7^l0g((A)" + (n2-l)(A) 

1 — a K^n-^' n^' 



^ log(n2)i-" 



1 — a 
= 21ogn. (25) 

Now we notice that the Renyi entropy is independent of a and furthermore coincides with the von 
Neumann entropy. This means that the reduced density matrix of a large block is proportional 
to n^-dimensional identity matrix. In other words, a sufficiently large block of neighbouring spins 
in our SU(n) VBS ground state is maximally entangled with the rest of the chain. Finally, we 
consider the analytic property of the Renyi entropy on the complex-a plane. In the large block 
limit, Sa{poo) is independent of a and hence Sa is completely analytic on a-plane. On the other 
hand, if we consider the finite-size block, Sa{pL) has branch cuts starting from Obc- The branch 
points Obc are determined from the conditionAo,o(-^)" + ("-^ — l)Ai,m7^o,o(-^)" = 0- The explicit 
value of Obc is given by 

^ ±(2m + l)7ri + log(n^ - 1) 
""^ log(Ao,o(i)/A;,„^^o,o(i)) ^ ' 

with m G Z. Since all the eigenvalues Xi^m become degenerate in the limit of large block sizes, 
the denominator of Eq. (l26p becomes zero and hence Obc converge to the point at infinity. This is 
the reason why the Renyi entropy in this limit is completely analytic on the complex a-plane. We 
also find the novel even-odd alternation of the real part of Obc, ^-e., Re[abc] > for even L, while 
Re[abc] < for odd L. 

C. Reduced density matrix as a projector onto the subspace of edge states 

As we have seen in the previous subsection, the von Neumann and Renyi entropies are obtained 
from the reduced density matrix for end two spins {pij- The obtained results indicate that the 



10 



!<!'„ 




11 2 2 



L L 



FIG. 2: Construction of the edge state |p, q) in an open boundary spin chain. A white (black) dot repre- 
sents the SU(n) fundamental (conjugate) representation. A circle denotes the projection onto the adjoint 
representation. The dotted line corresponds to the state \^p,q)- 



block of bulk spins is maximally entangled with the rest in the limit of large block sizes. The 
number of degrees of freedom in the subsystem can be counted from the von Neumann entropy as 
n?. This number coincides with the number of edge states which are degenerate ground states of 
the AKLT model with an open boundary condition. In the case of SU(2) AKLT model, the close 
relation between the von Neumann entropy and the number of edge states has been extensively 
discussed 



35 



46[. In this subsection, we elucidate the direct relation between the reduced density 
matrix ()14p and the edge states. First, we consider the open boundary SU(n) AKLT model with 
L sites. The Hamiltonian is given by -ffopen = X]fc=i ^k- The only difference from Eq.(l6|) is that 
there are no boundary terms. The basis of the edge states is constructed as: 



\P^(l) 



a 



p,<i 



E 

(ti.mi) 



Yl \^h,-rm)-U---\'^l 



L-i-mL-il L-lL-l 



#(0,0) 



X PLLiiUp,qUl,,m,...Ui 



L-i,nii_i 



^)l^o,o) 



LLJ^ 



(27) 



where Cp^q is a normalization factor. Any linear combination of (|27p is apparently the ground state 
of //open- The graphical representation of the construction of this state is shown in Fig. ([2]). The 



following orthogonality relation of edge state holds: {p, q\r, s) = Cpg{n'^ 



l)5p^rSq,s^-p~qiL) (the 



11 



subscripts —p and —q are modulo n). This can be shown as follows: 
{p, q\r, s) = Cp^qCr,s ^ ■■■ X] 

7^(0,0) 7^(0,0) 

7^(0,0) 7^(0,0) 7^(0,0) 

X LL{^l,m\{Ur,s ^ I)(.Ul,,n,,...Ul^_^,n,^_^ ^ I)\^0,0)lL 

= Cp^qCr,s{n^ - 1)^"^ 2Z '^l\rn'{L - l)i£($p+/'.q+m' | (1 - |$o,o)LL(^0,o|)|^r+/',s+m')LL 

{V,m') 

= C72g(n2-l)^,5p,,5g,,A„p,_,(L). (28) 

Here we have recalled Eq. (116p and have used the relation 

(!l,mi) (!£_l,™£_l) 

7^(0,0) 7^(0,0) 

= (n2 - i)^-Vrrr = ("' - 1)''"' E ^'.™(^ - i)l^^m)LL(^^"^l- (29) 

(/,m) 



The explicit form of the normalization factors Cp^q are given by Cp^q = Yj \J {r? — X)^\-p-q(V). 
Next, we try to write pi^ in terms of the basis of edge states. By the original definition, 

Pi^ = Tro,iTTlG')(G| 

~ C^2 _ i\L 2^ 2^ 2-^ ■■■ 2^ 2^ 

^(O'O) 7^(0,0) #(0,0) ^(0,0) 

X \'^h,-m,){'^l[,-m[\ ■ ■ ■ \^lL-u-mL-i){^l'^_„-m'^_,\ 

Then comparing with Eq. ()27p . we obtain 

(P.9) 

where A_p__g(= A„_p,i_g) was defined in Eq. (p^ . Therefore we can conclude that the reduced 
density matrix of a block of contiguous spins in the ground state is completely characterized by 
the ground states of corresponding open spin chain. In the limit of large block sizes, i.e., L -^ oo, 
PL can be written as 

In this limit, the limiting density matrix can be regarded as a projector which projects on a subspace 
spanned by a degenerate ground states for open boundary AKLT model. 



12 



III. PERIODIC SU(n) VBS STATE 

In this section, we shall focus on the von Neumann and the Renyi entropies of the periodic SU(n) 
VBS state. The periodic VBS state with length N can be constructed by acting the projection 
operator on the edges of the VBS state with boundary spins. We use the following representation 
for the ground state. 

10.) = ^ E - E E ■■■ E 

7^(0,0) ^(0,0) 7^(0,0) 5^(0,0) 

X \^h,-mi)ll ■ ■ ■ l^^ii^i.-mi^i/X^L-ll^'i+i.-mL+i/L+TL+l ' ' " l^'iV,-m]v)lVAf 

X ^LL(f^«L+l,mL+i---f^U,mjv^ii,mi---f^«L-i,mL-l®-^)I^O,0>LX, (33) 

where M"^ = {n? — 1)^(1 + {n^ — l)pn{N))/v? is a normalization factor. In the above state, the 
L-th site consists of the left (D) and right (D) "spins" in the original open SU(n) VBS state. The 
reduced density matrix of a block of contiguous spins with length L is given by 

Pn,l = Tri+i...7v|Gp)(Gp| 

~ j^ 2^ 2^ '" 2^ 2^ 2^ " ' 2^ 

(il.mi) ((' m') (lL-l,mx^_i) (l' m' ) (Il + 1-"^L + i) ('jV-^Jv) 

7^(0>0) ^(0,0) #(0,0) "- _i(o,o) #(0,0) #(0,0) 

X \'^h-mi)ll{^l[-m[\ ■ ■ ■ \^lL-i,-mL-i)L:^L-li^l'L^j-m'^_-^\ 

X P^,x(^i^3 ® I)\^ofl)LL{'^o,o\iV(Us /)tp^^, (34) 

where Vi = Ui^^rm ■ ■ ■ C^«i,_i,m£_i, Vl = Ui'^^^n', ■ ■ ■ Uv^_^,m'^_^ and U^, = Ui^^^^rnL+i ■ ■ ■ Uij^,^^. Here 
it is convenient to introduce the following density matrix: 

Ol+1'"^l+i) {'jv.»"jv) 

#(0,0) #(0,0) 

= Y. ^l^rniN - L)\<^l,m)LL{^l,m\. (35) 

{l,m.) 

The above expression is obtained from Eq. (J16p by replacing L with N — L. Eq. (|34p can be written 
in terms of Pjrzr as: 

PN,L = -^(n^-l)^-^ E E ••• E E \'^h,-m,)lli'^l'„-m[\ 

T^C'O) #(0,0) #(0,0) ^(0,0) 

= ^in'-l)''-'T.^UN-^) E E ••• E E l^h,-™.>n(^/i,-sl 

(/,m) (il,™l) (l{,m[) (lL-l,mL_i) (!^_-^ ,m'^_^) 

7^(0.0) ^(0,0) #(0,0) ^(o_o) 

X ••• |$^,_„-„^,_,)xrTi_l($^i_,,-™',_J^L^(^l®^)l'^/,™>LL(^Mnl(Vl'®/)V^^. (36) 
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Here we recall the basis of edge states (j27|l and obtain 

PN,L = jr^in -1) 2^ -^^2 \l,m){l,m\, 

(l,m) ''"» 

= Y.Xi,m{N,L)\l,m){l,ml (37) 

{l,m) 

where Xi^rn{N, L) = ''^2r<'i — ^(^^ ~ i^n-l ^^^^ eigenvalues of the reduced density matrix. They are 



explicitly given by 



J^ (l+(n^-lK(A.L))(l+(n-l)p4^)) (Z,^) = (0,0) 



Putting n = 2 into Eq. (f38l) reproduces the results of 5" = 1 case obtained in Ref.[[4a]]. 

Now we know the reduced density matrix of a block with arbitrary system size A^, subsystem size 
L and internal degrees of freedom n, we can study the finite size effect on entanglement properties. 
In the limit of A^ ^ cxd, the set of eigenvalues of the density matrix pn,l becomes equivalent to 
that of the VBS state with boundary spins, and hence the entanglement properties are identical 
to each other. The explicit form of the entanglement is similar to the one in previous section. The 
von Neumann and the Renyi entropies for periodic case are explicitly written as 

SipN,L) = -Ao,o(iV, L) log Ao,o(iV, L) - (n^ - l)Az,^^o,o(A', L) log Xi,m^o,o{N, L), 
Sa{pN,L) = -^log{Xo,o{N,L) + {n^-l)Xi,m^ofi{N,L)). (39) 

Similarly to the previous section, both S{pn,l) and Sa{pN,L) approach 2 log n in the thermodynamic 
limit. Here what we mean by thermodynamic limit is the limit in which A^ ^ oo, L — > oo with 
parameter L/N be fixed. 

As we saw in the calculation of the reduced density matrix, we found non-trivial consequence 
that the reduced density matrix is completely written in terms of the edge states of the open 
boundary AKLT model which corresponds to the Hamiltonian describing the subsystem ripped off 



by tracing out. It clarifies the ed, 
discussed in several papers 



35 



^e state interpretation of the entanglement entropy, which is also 



45 



46( 1 ■ in more detail. We can see the separation of the eigenvalues 
of the reduced density matrix, i.e., Xofi{N,L) ^ Xi^rn^ofi{N,L) in finite size systems. Here we shall 
explain the physical meaning of this separation. This is due to qualitative difference between the 
singlet state |0,0) and the adjoint states |/,m), (/,m) ^ (0,0) induced by the effective coupling of 
residual edge "spins" on the boundaries of the subsystem. In the thermodynamic limit, the reduced 
density matrix is proportional to the identity matrix, which indicates the absence of the effective 
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coupling in this limit and hence the edge spins behave freely without interaction. Therefore, the 
von Neumann entropy is proportional to the logarithm of the number of degerees of freedom due 
to the residual edge spins. 

IV. SUMMARY AND DISCUSSIONS 

We analyzed entanglement in the ground state of the SU(n) version of AKLT model. We 
consider a block of spins in the ground state, it is in the mixed state. We evaluated the von 
Neumann entropy and the Renyi entropy of the block. We first examined the VBS ground state 
with boundary spins. We found that the great simplification occurs in the limit of large block sizes. 
In this case the Renyi entropy is independent of the parameter a and furthermore it coincides with 
the value of von Neumann entropy 2 Inn. This means that the density matrix of the block is 
proportional to identical matrix of dimension n^. We clarified that subspace of eigenvectors of 
the density matrix with non-zero eigenvalues describes the degenerate ground state of the block, 
i.e., edge states. We studied the finite size effect on the Renyi entropy in terms of analyticity on 
the complex a-plane. Then we studied the periodic VBS state. In this case, we also obtained 
the exact expression for the reduced density matrix and found the essential simplification in the 
thermodynamic limit. 

V. ACKNOWLEDGEMENTS 

The authors are grateful to S. Murakami, and Y. Hatsugai for fruitful discussions. This work was 
supported in part by NSF grant DMS-0503712 (V.E.K.), Grant-in- Aids (Grant No. 15104006, No. 
16076205, and No. 17105002) and NAREGI Nanoscience Project from the Ministry of Education, 
Culture, Sports, Science, and Technology. H.K. was supported by the Japan Society for the 
Promotion of Science. 



[1] A. Osterloh, L. Amico, G. Falci, and R. Fazio, Nature (London) 416, 608 (2002) 

|arXiv:quant-ph/0202029| . 
[2] T. J. Osborne, M. A. Nielsen, Phys. Rev. A 66,032110(2002) ' arXiv:quant-ph/020216^ . 
[3] G. Vidal, J. I. Latorre, E. Rico and A. Kitaev Phys. Rev. Lett. 90, 227902 (2003) 

|arXiv:quant-ph/0211074] . 
[4] J. I. Latorre, E. Rico, and G. Vidal, Quant. Inf. Comput. 4, 048 (2004) [arXiv:quant-ph/030409^ . 



15 

[5] B.-Q. Jin, V. E. Korcpin, J. Stat. Phys. 116, 79 (2004) |arXiv:quant-ph/0304108| . 

[6] A. R. Its, B.-Q. Jin, V. E. Korcpin, J. Phys. A: Math. Gen. 38 2975 (2005) [arXiv:quant-ph/0409027| . 

[7] M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 87, 017901 (2001) [arXiv:quant-ph/0009060] . 

[8] F. Verstraete, M. A. Martin-Dclgado, and J. I. Cirac, Phys. Rev. Lett. 92, 087201 (2004) 
|arXiv:quant-ph/031108';1 . 

[9] J. K. Pachos, M. B. Plcnio, Phys. Rev. Lett. 93, 056402 (2004) |arXiv:quant-ph/040110^ . 

[10] H. Fan and S. Lloyd, J. Phys. A 38, 5285 (2005) |arXiv:quant-ph/0405130] . 

[11] V. Vedral, New J. Phys. 6 102 (2004), [arXiv:quant-ph/0405102l . 

[12] Y. Chen, P. Zanardi, Z. D. Wang, and F. C. Zhang, New J. Phys 8, 97 (2006) [arXiv:quant-ph/040722^ . 

[13] A. Hamma, R. lonicioiu, and P. Zanardi, Phys. Lett. A 337, 22 (2005) |arXiv:quant-ph/0406202l . 

[14] V. Popkov and M. Salerno, Phys. Rev. A 71, 012301 (2005) [arXiv:quant-ph/040402^ 

[15] J. P. Keating and F. Mezzadri, Commun. Math. Phys, 252, 543 (2004) |arXiv:quant-ph/0407047| . 

[16] S-J. Gu, S-S. Deng, Y-Q. Li, and H-Q. Lin, Phys. Rev. Lett. 93, 086402 (2004) 



|arXiv:quant-ph/0405067i 
[17] S-J. Gu, H-Q. Lin, and Y-Q. Li, Phys. Rev. A 68, 042330 (2003) |arXiv:quant-ph/0307131| . 
[18] J. Wang and S. Kais, Phys. Rev. A 70, 022301 (2004) [arXiv:quant-ph/0405087| . 
[19] V.E.Korepin, Phys.Rev. Lett. 92, 096402 (2004) [arXiv:cond-niat/031105"^ . 
[20] K. Audenaert, J. Eisert, M. B. Plenio, and R. F. Werner, Phys. Rev. A 66, 042327 (2002) 



|arXiv:quant-ph/0205025| . 
[21] C. H. Bennett, D. P. DiVincenzo, Nature 404, 247 (2000) . 
[22] S. Lloyd, Science 261,1569(1993); ibid 263, 695(1994) . 

[23] M. Levin and X. G. Wen, Phys. Rev. Lett. 96, 110405 (2006) jarXiv:cond-mat/051061^ . 
[24] A. Kitaev and J. PreskiU, Phys. Rev. Lett. 96, 110404 (2006) [arXiv:hep-th/0510092] . 



[25] A. Renyi, Probability Theory, (North-Holland, Amsterdam, 1970). 

[26] A.Afflcck, T.Kennedy, E.H.Licb and H.Tasaki, Commun. Math. Phys. 115, 477 (1988) 

[27] A.Affleck, T.Kennedy, E.H.Licb and H.Tasaki, Phys. Rev. Lett. 59, 799 (1987). 

[28] F.D.M.Haldane, Phys. Lett. 93A, 464 (1983); Phys. Rev. Lett. 50, 1153 (1983). 

[29] A.N.KiriUov and V.E.Korepin, ALGEBRA AND ANALYSIS [St. Petersburg Mathematical Journal], 

1,47 (1990). 

[30] J. J. Garcia-Ripoll, M. A. Marti'n-Delgado, and J. I. Cirac, Phys. Rev. Lett. 93, 250405 (2004) 



|arXiv:cond-mat/0404566| . 
[31] F. Verstraete and J. I. Cirac, Phys. Rev. A 70, 060302(R) (2004) |arXiv:quant-ph/0311130| . 
[32] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

[33] D. P. Arovas, A. Auerbach, and F. D. M. Haldane, Phys. Rev. Lett. 60, 531 (1988). 
[34] H. Fan, V. Korepin, and V. Roychowdhury, Phys. Rev. Lett. 93, 227203 (2004) 



|arXiv:quant-ph/0406067| . 
[35] H. Katsura, T. Hirano and Y. Hatsugai, Phys. Rev. B 76, 012401 (2007) |arXiv:cond-mat/0702196| . 



16 

[36] M. Greiter, S. Rachel, and D. Schuricht, Phys. Rev. 75, 060401(R) (2007) [arXiv:arXiv:cond- 

mat/0701354]. 
[37] M. Greiter and S. Rachel, Phys. Rev. B 75, 184441 (2007) |arXiv:cond-mat/0702443| . 
[38] A. R. Its, B.-Q. Jin, and V. E. Korepin, [arXiv:quant-ph/060617^ . 
[39] F. Franchini, A. R. Its, B.-Q. Jin, and V. E. Korepin, [arXiv:quant-ph/060624^ . 
[40] F. Franchini, A. R. Its, B.-Q. Jin, and V. E. Korepin, J. Phys. A: Math. Theor. 40, 8467 (2007) 



|arXiv:quant-ph/0609098| . 

[41] F. Franchini, A. R. Its, and V. E. Korepin, |arXiv:0707!253^ . 

[42] Entanglement of SU(n) singlet formed by n fundamental representations is considered in the following 

article: Z. Nussinov and G. Ortiz, [arXiv:cond-mat/0702377] . 
[43] H. Fan, Phys. Rev. Lett, 92, 177905 (2004) |arXiv:quant-ph/031102^ . 
[44] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, (Cambridge Univ. 

Press, Cambridge, 2000). 
[45] S. Ryu and Y. Hatsugai, Phys. Rev. B 73, 245115 (2006) |arXiv:cond-mat/060123"^ . 
[46] T. Hirano and Y. Hatsugai, J. Phys. Soc. Jpn. 76, 074603 (2007) [arXiv:cond-mat/0703642] . 



